PASJ: Publ. Astron. Soc. Japan , 1-??, 
© 2012. Astronomical Society of Japan. 



Early evolution of spherical ejecta expanding into the circumstellar matter at 

ultra-relativistic speeds 

Toshikazu Shigeyama/ Akihiro Suzuki^'^, Ko Nakamura^ 
^Research Center for the Early Universe, Graduate School of Science, University of Tokyo, Bunkyo-ku, Tokyo 113-0033 
^Department of Astronomy, Graduate School of Science, University of Tokyo, Bunkyo-ku, Tokyo 113-0033, Japan 
^National Astronomical Observatory of Japan, Mitaka, Tokyo 181-8588, Japan 

(Received ; accepted ) 
Abstract 

We present a new self-similar solution describing early evolution of an ultra-relativistic flow resulting from a col- 
lision of homologously expanding spherical ejecta with the circumstellar matter, in which a shock wave propagates 
in the circumstellar matter while a weak discontinuity propagates in the ejecta at the sound speed. 
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1. Introduction 

Ultra-relativistic flows play important roles in some astro- 
physical phenomena like 7-ray bursts where fire balls with 
the Lorentz factors exceeding 100 are thought to drive the 
phenomena. There are several self-similar solutions involv- 
ing such ultra-relativistic flows. Blandford & McKee (1976) 
found self-similar solutions for flows resulting from point ex- 
plosions in the relativistic limit. Adopting their prescrip- 
tion, Nakayama & Shigeyama (2005) described the shock 
emergence from a plane-parallel atmosphere. Furthermore, 
Nakamura & Shigeyama (2006) found relativistic version of 
the spherical flow relevant to young supernova remnants in 
which ejecta interact with the circumstellar matter (CSM) 
through two shock waves and the contact surface in between. 
Here we present a new self-similar solution describing flows 
associated with ejecta expanding into the CSM in somewhat 
different manner. A shock wave propagating in the CSM leads 
the jet and forms a dense shell detaching itself from the homol- 
ogously expanding ejecta. The inner wave front is a rarefac- 
tion wave propagating into the ejecta. This phase follows the 
spreading phase where the dense shell expands almost freely 
(Meszaros & Rees 1993) and precedes the phase described by 
Nakamura & Shigeyama (2006). 

A similar flow has been observed in 2-D axisymmetric sim- 
ulations for jets driven by aspherical supernova explosions 
(Nagakura et al. 201 1; Suzuki & Shigeyama 201 1) though de- 
tailed analyses of the structure of the jet in the radial direction 
is prohibited by coarse grids to cover the large dynamic range. 
It is expected from the same calculation that the propagation 
of the ultra-relativistic jet in the vicinity of the axis of symme- 
try is well approximated by the corresponding flow assuming 
spherical symmetry because the fluxes of mass, momentum, 
and energy as well as the pressure gradient in the lateral di- 
rection are too small to affect the structure and propagation of 
the jet. Therefore a spherically symmetric solution presented 
in this paper will be useful to investigate the radial structure of 
jets near the axis of symmetry. 

In the next section, we discuss initial conditions to realize 
the flow described by the self-similar solution presented in this 



paper. In section 3, we present the basic equations governing 
the ultra-relativistic flows in spherical symmetry. In section 4, 
we describe the density and pressure distributions of the ejecta. 
In section 5, we formulate the flow propagating in the ejecta 
using a prescription similar to that presented in Nakamura & 
Shigeyama (2006). In section 6, we discuss the propagation of 
the shock in the CSM using the self-similar solution found by 
Blandford & McKee (1976). The boundary conditions at the 
contact surface are given in section 7. Section 8 presents the 
results. In section 9, we conclude the paper. 

2. Initial conditions 

In the non-relativistic limit, a collision of supersonic ejecta 
with the CSM results in two shock waves with a contact sur- 
face in between if the relative velocity Ve satisfies the inequality 
(Landau & Lifshitz 1987), 

\l 2pcSM 

where Pg denotes the pressure in the ejecta, pcsM the proper 
mass density of the CSM, and 7 the adiabatic index, which is 
assumed as 7 = 4/3 in the rest of the paper. This inequality usu- 
ally holds in supernovae. In reality, young supernova remnants 
like Cassiopeia A are known to have such a structure (Gotthelf 
et al. 2001). On the other hand, in the ultra-relativistic limit, a 
similar argument yields a different condition for leading to two 
shock waves, that is. 

Here the velocity is measured in units of the speed of light c. 

To validate this condition, we solve the shock tube problems 
with parameters satisfying and not satisfying the above condi- 
tion and present the results in Figure 1 . In the left panel, two 
fluids with the density p — 10'^ are in contact at the position 
a; = at the time t — 0. The fluid in the region a; < having 
the pressure P = 1 moves in the positive x direction with the 
Lorentz factor 7e = 100, while the fluid in the region a; > is at 
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Fig. 1. Solutions for a relativistic shock tube problem with two different sets of initial conditions. Distributions of pressure (dashed line), ve- 
locity (dash-dotted line), and density (solid line) at time t = 0.5 are shown. The initial conditions at t = are as follows. The density 
is equal to 10^. The fluid in the region with a; < moves with the Lorentz factor 100 in the positive x direction and has the pressure 
equal to 1. The pressure in x > is 10^ in the left panel while 2 X 10^ in the right panel. The pressure is scaled by a factor of lO"'^. 



rest and has the pressure P ~ 10^. Therefore these satisfy the 
above condition (2) for two shock waves and the resultant flow 
in the figure at i = 0.5 actually shows two shock waves. Here 
the plotted pressure is scaled by a factor of 10^^ and the veloc- 
ity is measured in units of the speed of light. If the pressure in 
the right fluid is enhanced by a factor of two, the initial config- 
uration does not satisfy the condition for two shock waves and 
results in the flow pattern presented in the right panel at < = 0.5. 
However, it should be noted here that the criterion (2) does not 
accurately define the boundary of these two flow patterns and 
should be regarded as an approximate guide line to guess the 
resulting flow pattern from initial conditions. 

The value of the right hand side of this condition (2) could 
be estimated for a collision of ejecta with a stationary stellar 
wind in spherical geometry as 
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where ri3 denotes the distance from the explosion site in units 
of 10^^ cm, T7 the temperature of the ejecta in units of 10^ K, 
the wind velocity in units of 10^ cm s^^ blown by the pro- 
genitor star, and AI5 is the mass loss rate in units of IO^^Mq 
yr^^. If this condition is relaxed due to the decrease of the 
temperature when the forward shock propagates outward, the 
reverse shock will be eventually generated. Therefore we will 
present a self-similar solution for the spherically symmetric 
flow composed of a shock wave propagating in the CSM and a 
rarefaction wave in the ejecta resulting from the collision. This 
phase precedes the reverse shock phase described by the so- 
lution of Nakamura & Shigeyama (2006). Since the scale of 
the structure is very small, i.e., on the order of r/F^ where F 
denotes the Lorentz factor of the shock wave, this solution de- 
scribes fine structures of axisymmetric jets in the radial direc- 
tion near the symmetry axis that multi-dimensional numerical 
simulations hardly resolve. 



3. Ultra-relativistic flow 

The spherically symmetric flow moving at speed v is gov- 
erned by the following equations in the relativistic limit 
(Blandford & McKee 1976), 



dt 

dln(p3^4) 



,dP 

Ik' 

4 dr^v 



dt r2 
dp' 1 dr^p'v 
dt dr 



dr 
= 0. 



(4) 



(5) 



(6) 



Here r and t denote the radial and time coordinates. The 
Lagrange derivative is described by d/dt. The pressure is 
denoted by P and the ultra-relativistic equation of state e = 
3P/p >> 1 has been assumed. The mass density in the fixed 
frame is denoted by p', which is related to the mass density p in 
the co-moving frame as p' = •jp, where 7 is the Lorentz factor 
of the fluid element. 

In constructing the self-similar solution, we follow the pro- 
cedure of Nakamura & Shigeyama (2006) with different inner 
boundary conditions. 

4. Homologous ejecta 

Homologously expanding ejecta resulting from a shock 
emerging from the stellar surface have a power law structure in 
terms of the Lorentz factor of the fluid element as Nakayama 
& Shigeyama (2005) and Kikuchi & Shigeyama (2007) clearly 
showed. The velocity of homologous ejecta at distance r 
and time t is defined as. 
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The corresponding Lorentz factor is 
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The mass density of the ejecta follows the power law in terms 
of the Lorentz factor as, 

n 

(9) 



where a non-dimensional parameter n and a constant a have 
been introduced. The pressure is assumed to have a polytropic 
expression with a varying coefficient K{r), 



■■K{r)ptl\r). 



(10) 



In Nakamura & Shigeyama (2006), the pressure in the homolo- 
gous ejecta was neglected because a strong shock wave propa- 
gates in the ejecta. In this paper, a rarefaction wave is assumed 
to propagate in the ejecta and we need to specify the pressure 
distribution. 

5. Rarefaction wave 

Suppose that the front of the rarefaction wave is located at 
r = i?i at time t. Then the boundary conditions are described 
as 



Pl7l =Pe(^l)7e(^l), 

P^=P, = K{R^)pt'\Rr). 
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The rarefaction front propagates at a speed Vi{t). Following 
the prescription of Blandford & McKee (1976), the corre- 
sponding Lorentz factor Fi is assumed to depend on time as 

Tl = At-"', (14) 

where A and m are constants. Therefore Ri is approximated 
by 

(15) 
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Using this expression, other variables in the ejecta at r = -Ri 
are expressed as 

7f(i?i) = (m + l)r?, 



Pe(-Rl 



(16) 

(17) 



The coefficient K{Ri) is assumed to be K{Ri) — KQ{t/toy 
where Kq, to, and / are constants. Using these di- 
mensional variables at the front, non-dimensional variables 
G{0, H{0 are introduced as follows. 

P{r,t) = K{R,)p*J'iRi)FiO, (18) 
72(r,t) = (m+l)r?G(0, (19) 
p{r, thir, t) = p,{RMRi)H{^). (20) 

Here a non-dimensional coordinate was defined as ^ = [1 + 
2(m -|- l)ri](l — r/t) so that the rarefaction front is located at 
^ = 1. Since the rarefaction wave propagates at the sound speed 
equal to 1 /\/3 relative to the matter in the relativistic limit, the 
Lorentz factors of the rarefaction front and the ejecta there are 
related with each other as 
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Thus the parameter m is evaluated as m = 1 /\/3- 

The conversion of coordinate system (r, t) to {£,, Fi) is de- 
scribed as 



*|:=-[l + 2(m+l)F?]|, 
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If we retain only leading terms, the resultant hydrodynamical 
equations deduced from equations (4)-(6) become three ordi- 
nary differential equations. 

dlnF 12[l-{m + l)^G]dlnG 



d^ 6[l + (m-|-l)CG] d^ 
[m(4n-12)+6?-24]G 
^ 6[l + (m+l)eG] ' 
8(to + l)(m + 2)^0 + 2m(4n - 6) 
~ ^{[l + {m + 1)^G]2 - 3[1 -{m + 1)CG]2} 

12; -4m -32 
4{[1 + (m + 1)CG]2 - 3[1 - (m + l)^G]^} ' 



dG 



dH 
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)dlnG 
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(25) 
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6. Flow in tlie sliocked circumstellar matter 

Finally, we will present equations governing the self-similar 
flow in the shocked CSM and the associated boundary con- 
ditions. These are identical to those in Blandford & McKee 
(1976). 

6.1. Shock conditions in the CSM 

The mass density, Lorentz factor, and pressure change across 
a strong shock wave propagating at a speed V by the following 
relations. 



P2 



E P272 = 2piF^ 
2 

72 = Y' 
n 2piF2 



(27) 
(28) 

(29) 



where the subscript 1 refers to the values in the shocked fluid 
at the shock front, the subscript i the values in the un-shocked 
external medium at the shock front. F = 1/ Vl — denotes 
the Lorentz factor of the shock front. To maintain the self- 
similarity, the Lorentz factor of the shock front F needs to 
evolve as Fi and is assumed to depend on time as F^ = Bt^™. 
Here B is a constant. The mass density in the CSM is assumed 
to have a power law distribution as p\ = br~'' where b and k 
are constants. Non-dimensional self-similar variables for the 
pressure /(x)> the Lorentz factor g{x)' and the mass density 
/i(x) are defined as 
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where the similarity variable x has been introduced as 
X={l + 2(m + l)r2}(l-0. 

6.2. Self-similar flow 

By substituting the expressions (30)-(32) into equations (4)- 
(6), equations governing the self-similar flow in the shocked 
CSM are obtained as 



din/ _ 4{2(m - 1) + k} - {m + k - 4)gx 



(m - 
(7m - 



gdx 

ding 

gdx (m + l)(4- 
dlnh _ 2(9m + 5fc-8) 

gdx 



-m-Sgx + 
3fc - 4) - (m 



■2)<?X 



85X + 5^X^) 
-2{5m + Ak-6)gx 



(m + l)(4-85X + ffV)(2-5X) 

{m + k-2)g^X^ 



(34) 
(35) 

(36) 



(m + 1)(4 - 8gx + 5'x2)(2 - ffx) ' 

Because m = 1/ V3, the similarity solution is uniquely deter- 
mined by specifying the parameter k. 

7. The flow at the contact surface 

At the contact surface, the density distribution has a discon- 
tinuity while the pressure and the velocity are continuously dis- 
tributed. Since the pressures in equations (30) and (13) need to 
evolve with the same power of time t, we obtain the relation 



Tn{2n + 3) + 3{k + l) = 12. 



(37) 



Using this relation in equations (24) and (25), the dependences 
of these equations on the parameters / and n are eliminated. 
The continuous distribution of pressure at the contact surface 
also relates the introduced constants as 



3F(Ce) 
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and 5f(Xc)Xc = 2, 



(39) 



The denominators of (36) and (26) indicate 
1 

TO + 1 

where and Xc denote the coordinates of the contact surface. 

From the continuity of the Lorentz factors at the contact sur- 
face, equations (3 1) and (25) yield the ratio of Lorentz factors 
of the two wave fronts as 



r 
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The asymptotic behavior of the density near the contact surface 

is described as 



/i«(Xc-x)"^'"+^-^"^, 

n 4(n + l) 

^0C(^ — ^c) 3%/3fc + l-12v^ 



(41) 
(42) 



The former equation indicates that the density in the CSM 
diverges if l/\/3 < < 4 — a/3/9, otherwise the density in 



the CSM becomes zero at the contact surface. From the lat- 
ter equation, when k = 0, i.e., the CSM has a uniform den- 
sity, the density in the ejecta diverges at the contact surface if 
n < 6\/3 - 3/2. When k = 2, i.e., the CSM is formed by a 
stationary stellar wind, the density in the ejecta diverges at the 
contact surface if n < 3a/3 — 3/2 . 

8. Results 

Results with a combination of parameters n = k = 2 are 
presented in Figure 2. The other parameters are as follows: 
a = 10*^, b = 5.02 X 10", Kq = 5.43 x lO^^^ and tn = 10^ 
in cgs units. These values are obtained from ejecta having the 
Lorentz factor of 300, the temperature of 10^ K, and the den- 
sity of 10~^° g cm~^ at the front. The ejecta collide atr = cto 
with the stationary wind with the mass loss rate of 10^ '' Mq 
yr~^. The wind velocity is assumed to be 1 x 10^ cm s~^. The 
plotted pressure, Lorentz factor, and density are normaUzed by 
those values at the shock front in the CSM. Note that the den- 
sity profile on the left side of the contact surface in this plot 
changes with time while the other profiles remain unchanged 
with time. 

From the obtained solutions, we can derive a relation be- 
tween the Lorentz factor and pressure of the ejecta and the den- 
sity of the CSM near the wave as 



7e 



/ 3F{^,)g{Xc}Po 
4/(Xc)G(ec)pcsM^ 



(43) 



by using equations (38) and (39). The factor 
F{QgiXc)/if{Xc)G{^c)) ranges from 2.2 to 4.5 for the 
parameter range of < fc < 2 with a fixed n{= 2), for example. 
Therefore the solutions satisfy the criterion (2) though we have 
rarefaction waves rather than reverse shock waves. Again, 
we should recognize the inequality (2) as an approximate 
guide Une to classify the flow patterns resulting from initial 
discontinuities partly because of the inhomogeneous structure 
of the flow. 

9. Conclusions and Discussion 

We have presented a new self- similar solution describing the 
early phase of collision between spherical ejecta expanding at 
ultra-relativistic speeds and the CSM. The flow is character- 
ized by the forward shock propagating in the CSM, the con- 
tact surface, and the rarefaction front with a weak discontinu- 
ity propagating in the ejecta. As the forward shock proceeds 
in the CSM, a reverse shock eventually forms. Then another 
self-similar solution by Nakamura & Shigeyama (2006) appro- 
priately describes the flow afterwards. Though these solutions 
assume spherical symmetry, they can be used to describe flows 
of relativistic jets in the vicinity of the axis of symmetry since 
the 2-D numerical simulations show that the flux in the trans- 
verse direction is negligible there. 

When the rarefaction described by the present solution is 
transformed into the reverse shock depends on the structures 
of ejecta and the CSM. If the ejecta are isothermal, the condi- 
tion (2) together with equation (3) suggests that the transition 
occurs at time t ^ 2 x 10^ sec. This corresponds to -^0.2 sec in 
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Fig. 2. Distributions of Lorentz factor (dasli-dotted line), density (solid 
line), and pressure (dashed line) normalized with those values at 
the shock front as functions of scaled distance from the contact 
surface. The parameters are n = 2, k = 2, 1 = 2 — 7\/3/9. 



the observer's frame. In reality, higher temperature or pressure 
in the inner ejecta will prolong the duration of this phase. In 
the self-similar solution presented in this paper, the distribution 
of pressure P{r, t) in the ejecta has a form of 

Poc 72-4^3^-4. (44) 

If the actual pressure distribution deviates from this and be- 
comes shallower, then the transition is expected to occur. The 
transition indicates a change of the time dependence of the 
shock Lorentz factor. It will be intriguing to investigate influ- 
ences of this change on the emission in connection with light 
curve features of some gamma-ray bursts. For example, GRB 
1 10205 A exhibits a steep rise in its optical light curve starting 
at a few hundred seconds after the detection while the X-ray 
light curve shows a small bump around the same period (Zheng 
et al. 20 11). This might be related to the transition taking place 
at r ^ 10^* cm from the explosion site and associated decel- 
eration of the shock front, though calculations of the emission 
from such a shocked relativistic matter is beyond the scope of 
the present paper. 
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